Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 5 February 2008 



(MN M?eX style file v2.2) 



Roche tomography of cataclysmic variables — III. Starspots 
on AE Aqr 

C. A. Watson/* V. S. Dhillon/ and T. Shahbaz^ 

^ Department of Physics and Astronomy, University of Sheffield, Sheffield S3 7RH, UK 
^ Instituto de Astrofisica de Canarias, 38200 La Laguna, Tenerife, Spain 



Submitted for publication in the Monthly Notices of the Royal Astronomical Society 

5 February 2008 

ABSTRACT 

We present a Roche tomography reconstruction of the secondary star in the cata- 
clysmic variable AE Aqr. The tomogram reveals several surface inhomogeneities that 
are due to the presence of large, cool starspots. In addition to a number of lower- 
latitude spots, the maps also show the presence of a large high latitude spot similar to 
that seen in Doppler images of rapidly-rotating isolated stars, and a relative paucity of 
spots at a latitude of 40°. In total, we estimate that some 18 per cent of the Northern 
hemisphere of AE Aqr is spotted. 

We have also applied the entropy landscape technique to determine accurate pa- 
rameters for the binary system. We obtain optimal masses of Mi = 0.74 Mq, M2 = 
0.50 M0, a systemic velocity 7 = -63 km s~^ and an orbital inclination of i = 66°. 

Given that this is the first study to successfully image starspots on the secondary 
star in a cataclysmic variable, we discuss the role that further studies of this kind may 
play in our understanding of these binaries. 

Key words: stars: novae, cataclysmic variables - stars: spots - stars: late-type - 
stars: imaging - stars: individual: AE Aqr - techniques: spectroscopic 



1 INTRODUCTION 

Cataclysmic variables (CVs) are short period binary sys- 
tems in which a (typically) late main-sequence star (the 
secondary) transfers material via Roche-lobe overflow to a 
white dwarf primary star. Although CVs are largely ob- 
served to study the fundamental astrophysical process of 
accretion, it is the Roche-lobe filling secondary stars them- 
selves that are key to our understanding of the origin, evo- 
lution and behaviour of this class of interacting binary. 

For instance, magnetic activity of the secondary star 
plays a pivotal role in the canonical theory of CV evolution 
by draining angular momentum from t he binary through a 
process known as magnetic braking (e.s. lKraf3ll967l iMestell 
ll968l.lSpruit fc Rittel^l983^ . This is beUeved to be responsi- 
ble for sustaining mass transfer between the components and 
to cause CVs to evolve to shorter orbital periods. At a sec- 
ondary star mass of ~ 0.25 Mq, corresponding to an orbital 
period of around 3 hours, the secondary star becomes fully 
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convective. This is thought to quench the dynamo mech- 
anism and shut down the mass transfer. Contact is later 
re-established at a shorter orbital period (~ 2 hours) due to 
angular momentum loss by gravitational radiation - explain- 
ing the dearth of CVs seen between 2-3 hou rs (the so-called 
period gap - see e.g. iTaam fc Spruidll989l). Although this 
is the canonical picture of CV evolution, lAndronov et alj 
(12003?) have noted that it is unclear if magnetic braking is 
the dominant angular momentum loss above the period gap, 
and there is no observational support for an abrupt change 
in angular momentum losses at the point where stars become 
fully convective. 

Furthermore, magnetic activity cycles on the secondary 
star have been invoked to ex plain the orbital period varia- 
tions observed for some CVs ilApplegatelll992l) . In addition, 
an increase in the number of magnetic flux tubes on the 
secondary star may cause the star to expand relative to its 
Roche lobe (e.g. lRichman et al.ll994|) and result in enhanced 
mass transfer - giving rise to an increased mass transfer rate 
through the accretion disc and a corresponding increase in 
the system luminosity. The additional mass flow would also 
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reduce the time required to build up enough material in the 
disc to trigger an outburst - resulting in shorter time inter- 
vals between consecutive outbursts. 

It is also speculated that the secondary star magnetic 
field could reach across to the disc and remove angular mo- 
mentum Irom its outer regions, causing material to accumu- 
late in the inner disc. An increase in the number of magnetic 
flux tubes on the secondary could then ca use m aterial to ac- 
cumulate in the inner regions faster (e.g. lAk et al. 20011 - 
resulting in shorter outburst durations. Thus, magnetic phe- 
nomena on the secondary star have been invoked to explain 
variations in the orbital periods, mean brightnesses, mean 
outburst frequencies, mean outburst durations and outburst 
shapes of CVs, as well as playing a major role in their evo- 
lution. 

From another (perhaps more important) perspective, 
by studying magnetic activity in CVs we can provide tests 
of stellar dynamo theories. Although faint (ttIu > 11), CVs 
contain some of the most rapidly rotating stars known, and 
therefore provide an excellent test-bed in which to study 
such models since they should have strong dynamo action 
(e.g. ^Rutten 1987). Furthermore, tidal forces are thought to 
have a strong i mpact on stellar dyn amos, suppressing differ- 
ential rotation iScharlemannl ll 982fl and forcing starspots t o 
form at preferred longitudes jHolzwarth fc Schussleill2003l) . 
Currently, observational studies of the impact of tidal effects 
on stellar dynamos are still lacking. Since the secondary stars 
in CVs are heavily distorted, studies of magnetic activity on 
these stars are highly desirable. 

Despite the obvious importance of magnetic activity to 
our understanding of the evolution and behaviour of these 
binaries, and the role that studying such activity on the 
secondary stars in CVs can play in furthering our under- 
standing of the solar-stellar connection, little work has been 
carried o ut in this field (w ith the notable exception of the 
work by IWebb et al.1 12002H . We aim to address this issue. 
In this paper we present the first images of starspots on 
the second ary star in a CV (AE Aqr) using Roche tomog- 
raphy (see 'Watson fc Dhillonll200lt iDhillon fc WatsonlbOOll 
and ■ Watson et aL.200 ^ for details') . 



2 OBSERVATIONS AND REDUCTION 

Simultaneous spectroscopic and photometric data of AE Aqr 
were acquired over two nights on 2001 August 9 and 10. The 
spectroscopic data were obtained using the 4.2-m William 
Herschel Telescope (WHT) and the simultaneous photome- 
try was carried out using the 1.0-m Jacobus Kapteyn Tele- 
scope (JKT). Both telescopes are part of the Isaac Newton 
Group of telescopes situated on the island of La Palma. 

2.1 Spectroscopy 

The spectroscopic observations of AE Aqr were carried out 
using the Utrecht Echelle Spectrograph (UES) and E31 grat- 
ing on the WHT. The SITel CCD chip with 2088 x 2120 
pixels (each of 24/im) was used. Centred at a wavelength 
of 5000A in order 114, this allowed orders 77-141 to be 
covered (corresponding to a wavelength range of 400lA - 
7475A, with significant wavelength overlap between adjacent 
orders). With a slit width of 1.11 arcsec (projecting onto 2.0 



pixels at the detector) a spectral resolution of around 46000 
(i.e. ~6.6 km s~^) was obtained. The majority of spectra 
were taken using a 200-s exposure time in order to minimise 
velocity smearing of the data due to the orbital motion of the 
secondary star. Comparison ThAr arc lamp exposures were 
taken approximately every hour to calibrate instrumental 
flexure. 

With this setup we obtained 88 usable spectra (with 
one 400-s exposure) on 2001 Aug 9, and 95 usable spec- 
tra on the following night. Both nights were clear (though 
dusty) with seeing varying between 1 and 1.4 arcsec at low 
air mass, worsening to 2 arcsec towards the end of the sec- 
ond night. The signal-to-noise of each spectrum peaked at 
~20-44 per pixel around 5580A. Table [H gives a journal of 
the observations undertaken at the WHT. 

2.1.1 Data reduction 

The raw data were first bias subtracted and the overscan 
regions cropped. Pixel-to-pixel variations were corrected for 
using a tungst en lamp flat -field. An order-by-order opti- 
mal extraction llHornelll986l) was then carried out using the 
PAMELA package written by torn Marsh. The arc spectra 
were extracted from the same locations on the detector as 
the target and fitted with a fourth-order polynomial, which 
gave a typical rms scatter of better than O.OOlA. After ex- 
traction and wavelength calibration, each order was then 
inspected visually for any obvious wavelength calibration 
errors and to ensure that no jumps were present between 
adjacent orders. Orders that were affected by bad columns 
on the SITel CCD chip were then discarded (these were or- 
ders 113 and 120 centred on 5055A and 4757A respectively). 

We should note here that Roche tomography cannot 
be performed on data which has not been corrected for slit 
losses. This is because the variable contribution of the sec- 
ondary star to the total light of a CV forces one to use rel- 
ative line fluxes during the mapping process and prohibits 
the usual method of normalising the spectra. We corrected 
for slit losses by dividing each AE Aqr spectrum by the ra- 
tio of the fiux in the spectrum (after folding the spectrum 
through the photometric filter response function) to the cor- 
responding photometric flux (see Section. 12.21 . 

2.2 Photometry 

Simultaneous photometric observations were carried out on 
the JKT using a Kitt Peak V-band flher and the SITe2 CCD 
chip with 2088x2120 pixels, each of 24/im. In order to re- 
duce readout time, 2 windows were created. The first window 
(measuring 920x1000 unbinned pixels) was setup to cover 
the target and the b r ightes t photometric standards from 
iHenden fc Honevcutd il995l) . and a second window (mea- 
suring 920x80 unbinned pixels) covered the bias strip. The 
chip was then binned by a factor of 2x2. 

2.2.1 Data reduction 

Since the master bias-frame showed that large scale struc- 
ture in the bias was insignificant, the bias level of the target 
frames was removed by subtracting the median of the pixels 
in the overscan region. Pixel-to-pixel variations were then 



Roche tomography of cataclysmic variables - III. Starspots on AE Aqr 3 



Table 1. Log of the WHT spectroscopic observations of AE Aqr, the spectral-type template stars, a telluric-correction star and the 
spectrophotometric standard HR 26. The first column gives the object name, columns 2-4 list the UT Date and the exposure start and 
end times, respectively, columns 5—6 list the exposure times and number of spectra taken for each object and the final column indicates 
the type of science frame taken. 



Object 


UT Date 




UT Start 


UT End 




No. spectra 


Comments 


AE Aqr 


2001 Aug 


09 


21:01 


21:55 


200 


11 


Target spectra 


AE Aqr 






21:57 


22:04 


400 


1 




AE Aqr 






22:05 


04:22 


200 


76 




HR 26 


2001 Aug 


10 


04:29 


04:29 


10 


1 


Telluric B9V star 


HR 26 






04:32 


04:35 


200 


1 




HR 26 






04:37 


04:42 


300 


1 




HD 19455 






04:56 


05:09 


300 


2 


Spectrophotometric standard 


HR 9038 






05:13 


05:25 


300 


2 


K3V template 


HR 8372 






05:37 


05:49 


300 


2 


K5V template 


HR 8881 






05:53 


06:05 


300 


2 


KIV template 


AE Aqr 






20:49 


04:37 


200 


95 


Target spectra 


HR 8086 


2001 Aug 


11 


04:44 


04:56 


300 


2 


K7V template 


HR 8382 






05:04 


05:16 


300 


2 


K2V template 


HR 8 






05:20 


05:31 


300 


2 


KOV template 


Gl 157A 






05:42 


06:14 


300 


4 


K4V Template 



removed by dividing the target frames through by a mas- 
ter twilight flat-field frame. Aperture photometry was per- 
formed using the package photom. Finally, the photometry 
data were corrected for atmospheric extinction by subtract- 
ing the magnitude of a nea rby comparison star (AE Aqr-2; 
iHenden fc Honevcuttll995h . The hght curves for both nights 
are shown in Fig. 

Inspection of the light-curves clearly shows the flaring 
activity so typical of AE Aqr (e.g. |BeskrmaiM^LM_a Jj|l99d; 
ISkidmore etal] 120031: iPearson. Home fc Skidmord I2003D . 
with flare amplitudes approaching ~0.5 mag. 



2.3 Continuum fitting 

As mentioned earlier in Section [2. 1.11 there is an unknown 
contribution to each spectrum from the accretion regions. 
The constantly changing continuum slope due to, for exam- 
ple, flaring (see Fig or from the varying aspect of the 
accretion regions means that the construction of, for in- 
stance, a master continuum flt to the data (such as that 
done bv lCollier Cameron fc Unruhlll99^ is not appropriate 
under these circumstances. Furthermore, since the contribu- 
tion of the secondary star to the total Ught of the system 
is constantly varying, normalisation of the continuum would 
also result in the photospheric absorption lines from the sec- 
ondary star varying in relative strength from one exposure 
to the next. Therefore, we are forced to subtract the contin- 
uum after each spectrum has been corrected for slit-losses 
using simultaneous photometry (Section 12.211 . 

The continuum fitting was carried out by placing spline 
knots at locations that were relatively line-free. The posi- 
tions of these spline knots were chosen on an individual 
spectrum-by-spectrum basis such that a smooth and visually 
acceptable fit was acquired, ignoring regions that included 
emission lines. Naturally, this is a very subjective procedure 
as it is particularly difficult to judge where the true contin- 
uum actually lies amongst the forest of heavily broadened 
lines. 



- D.2 0,3 0,4 0,5 0,6 0,7 0,6 0,9 

Orbital Phase 




2,6 2,7 2,8 2,9 3 3,1 3,2 3,3 



Orbitol Phase 

Figure 1. AE Aqr V-band JKT light curves. Top: data from 
2001 August 09. Bottom: data from 2001 August 10. The data 
have been phased according to the ephemeris derived in SectionlSl 
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In order to assess the impact of continuum mis-fits on 
our analysis we experimented with different fits, including 
those that were quite obviously incorrect. We were com- 
forted to find that a poor continuum fit did not adversely 
affect the shape of the line profiles after carrying out the 
least squares deconvolution (LSD) process described in Sec- 
tion |3 Evidently the results of this process are dominated 
by the sheer number of lines used rather than any contin- 
uum mis-fits, even if these mis-fits amount to quite major 
departures from the true continuum shape. 

We did find, however, that we systematically fit the con- 
tinuum at too high a level, leading to LSD profiles with con- 
tinuum (subtracted) regions lying significantly below zero. 
This was rather unexpected as the heavily broadened lines of 
the secondary act to 'depress' the observed continuum level, 
which would most likely lead to the continuum being fitted 
at too low a level. We can only surmise that we were suc- 
cessful in finding regions devoid of lines and that we fitted 
across the top of the noise in these regions. This was easily 
solved, however, by shifting the continuum fit to lower levels 
until the continuum in the LSD profiles lay at zero. This did 
not affect the shape of the line profiles obtained. 



3 LEAST SQUARES DECONVOLUTION 

Previous Roche tomograms have been restricted to using a 
single absorption or emission line on a 4-m class telescope 
coupled with an intermediate resolution spectrograph (e.g. 
LWat son ct al. 2003). Although successful at picking out the 
strong signatures of irradiation on the secondary stars, the 
features due to starspots are far more subtle. In order to 
achieve high enough signal-to-noise to detect starspot fea- 
tures in absorption-line profiles of the secondary in AE Aqr 
we have applied the technique of LSD. 

LSD effectively stacks the ~1000's of stellar absorp- 
tion lines observable in a single echelle spectrum to pro- 
duce a single 'average' profile of greatly increased signal- 
to-noise. Theoretically, the multiplex gain in signal-to-noise 
is the square root of the number of lines observed. The 
technique was first appl ied to spectropolari metric observa- 
tions of active stars by iDonati et alJ (Il997l) and has since 
been used in a large num ber of Doppler imaging studies 
by a variety of groups (e .g. 'Barnes ct al."l998'; 'Lister et^| 
19991: iBarnes et alJipOOol: iBarncs fc Collier Cameron, . 20011 
Barnes et al. 2001al:]Barnes et al. 2001bl|Jeffers et al.ll2002t 
Barnes et alJ l2004t iMarsden et alJ l2005l) . For supplemen- 
tary d etail s of the LSD te chnique, see also lCollier CameronI 
J200ll) and lBarned Jl999h . 

LSD assumes that starspots affect all of the rotation- 
ally broadened profiles in the same manner. Although the 
amplitude of the starspot bumps may vary from line to line, 
their morphology should be identical. LSD also requires the 
positions and relative strengths of the lines observed in each 
echelle spectrum to be known. Currently, we use line lists 
generated by the Vienna Atomic Line Database (VALD ) for 
this purpose (see iKupka et al]ll999l: iKupka et al]|200CD . 



The spectral type of AE Aqr has been determined 


to lie in the range K3-K5 V/IV JCrawford & Kraft 


1956; 


Tanzi. Chincarini & Tarenehll981ilWadell982i:iBruch 


,19911 


Welsh. Horne & Okd Il993l: iReinsch & BeuermannI 


1994 


Casares et al.lll996l). Therefore, a line-list for a stellar at- 



mosphere with Tefj — 4750K and logg = 4.5 (the closest 
approximation available in the database to a K4V spectral 
type) was downloaded for use in the LSD process. Since the 
line-list obtained from vald contain normalised line-depths, 
whereas Roche tomography uses continuum subtracted spec- 
tra, each line-depth was scaled by a fit to the continuum of 
a K4V template star. This meant that each line's relative 
depth was now correct for use with continuum subtracted 
data. 

Although A E A q r is slightly e volved 
jEracleous & Horn(^ Il996l: IWvnn. King fc Hornd Il997l : 

^carson ct al. 200^ and h as non-solar abundances due to 
possible CNO cyclin g (e.g. iMauche. Lee fc Kallnianl Il997l : 
ISchenker et al.ll2002h . our choice of line-lis t is unli k ely to 
affect the results presented in this paper. iBarnej lll999l) 

investigated the robustness of the LSD process with respect 
to the choice of the line-list spectral type used. By carrying 
out LSD on the rapidly rotating single star PZ Tel (spectral 
type KO V) usin g different line lists (from spectral type 
G2-K5), iBarnesI |599) found that the LSD process was 
insensitive to the use of a line-list of incorrect spectral type. 
He also concluded that, therefore, the effects of metallicity 
would also have little impact on the LSD process. 

One further complication in the LSD process that is 
possibly worth considering is whether changes in the con- 
tinuum lightcurve can result in significant changes in the 
tilt of the spectrum. This would lead to lines receiving a dif- 
ferent weighting at different orbital phases which may cause 
slight changes in the deconvolved profiles line-depths and 
shift the effective 'central wavelength' of the observations 
(which would alter the effective limb darkening coefficient 
seen). These effects are likely to be small, and certainly can- 
not be seen in this work. 

We carried out LSD over the wavelength range of 
4200A - 6820A. Bluer than ~4200A, the signal-to-noise is 
poor and degrades the quality of the final deconvolved pro- 
file. For regions redder than ~6820A, the number of pho- 
tospheric absorption lines is greatly diminished and telluric 
features become more abundant - inclusion of which would 
degrade the final deconvolved profile. We also masked out 
regions where strong emission lines were evident, such as 
around H7, H/3 and Ha. Altogether, a total of 3843 lines 
were used in the deconvolution process, resulting in an av- 
erage multiplex gain in signal-to-noise of 26 over single- 
line studies. Errors were propagated through the LSD pro- 
cess, which results in an under-estimate of the signal-to- 
noise rati o in t he de convol ved profiles (as a lso found by 
^arsdon et al.l ^00^ - see lWade eralll2000l for more de- 
tails). This results in our images being reconstructed to re- 
duced values less than 1, but does not affect the quality 
of the final maps. 

The individual deconvolved profiles obtained on 2001 
August 9 and 10 are displayed in Figs.|5|and ^ respectively. 
These clearly show bump features moving from blue (-ve 
velocities) to red (-|-ve velocities) across the profile. These 
features, as well as the orbital motion and variation in the 
projected equatorial rotation velocity vsini, are also obvi- 
ous when the LSD profiles are trailed (see Fig|lJ. In order 
to enhance the spot features, we have subtracted a theoret- 
ical profile from each LSD profile and removed the orbital 
motion - the result is shown in Fig. |^ (The theoretical pro- 
file removed from the data was generated using our Roche 
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velocity (km/s) 





-200 200 

velocity (l<m/s) 



Figure 4. A trail of the deconvolved profiles from observations on 
2001 August 9 (left) and 2001 August 10 (right). The small gaps 
in the phase coverage are at times when arc spectra were taken for 
the purpose of wavelength calibration. For clarity, gaps between 
individual exposures due to readout time are not shown. In these 
trailed spectra, features due to starspots will appear bright and 
several such features can be seen which appear on both nights 
(see Fig.|^for plots with increased contrast). Also evident is the 
variation in vsini showing the maximum at quadrature (phase 
0.25 and 0.75) due to the changing aspect of the tidally distorted 
star. 




J* J ■ . " .- 'A <i\7 'J'E' 



-200 -100 100 200 -200 -100 100 200 

velocity (km/s) velocity ^km/s.) 

Figure 5. A trail of the deconvolved profiles from observations 
on 2001 August 9 (left) and 2001 August 10 (right). The orbital 
motion has been removed assuming the binary parameters found 
in Sectionlsl which allows the individual starspot tracks across the 
profiles and the variation in v sin i to be more clearly observed. In 
order to increase the contrast and allow the starspot features in 
the trailed spectrum to be more readily seen, a theoretical profile 
(see text) has also been subtracted from each LSD profile before 
trailing. In this case, starspot features now appear dark in the 
plots. Once again, the small gaps in the phase coverage are at 
times when arc spectra for the purpose of wavelength calibration 
were taken and gaps between individual exposures due to readout 
time are not shown. 



tomography code adopting the binary parameters found for 
AE Aqr (see Section EJ, assuming a featureless stellar sur- 
face and the same limb-darkening coefficients as used in the 
final reconstructions). 



4 ROCHE TOMOGRAPHY 

Roche tomogr ap hy (|^u^«i^^Miillo^_^ 

DhiUon fc WatsonI 1200 it IWatson fc DhiUonI boOU 
W£feor^^dn20o3) is analogous to the Doppler imag- 
ing technique used to m ap rapidly rotating single stars 
(e.g. IVoet fc Penrodlll983l) . Specifically, Roche tomography 
is designed to map the Roche-lobe-filling donor stars in 
interacting binaries such as CVs and X-ray binaries. In 
Roche tomography, the secondary star is assumed to be 
Roche-lobe-filling, locked in synchronous rotation and to 
have a circularized orbit. The secondary star is then mod- 
elled as a series of quadrilateral tiles or surface elements, 
each of which is assigned a copy of the local (intrinsic) 
specific intensity profile. These profiles are then scaled to 
take into account the projected area of the surface element, 
limb darkening and obscuration, and then Doppler shifted 
according to the radial velocity of the surface element at 
that particular phase. Summing up the contributions from 
each element gives the rotationally broadened profile at 
that particular orbital phase. 

By iteratively varying the strengths of the profile con- 
tributed from each element, the 'inverse' of the above pro- 



cedure can be performed. The spectroscopic data is fit to 
a target s^nd a unique map is selected by employ- 
ing the maximum-entro py memsys algorithm developed by 
ISkilling fc Brvaiil (|l984h . Unless otherwise stated in the text, 
we employ a moving uniform default map, where each ele- 
ment in the default is set to the average value in the re- 
constructed map. Since the maximum-entropy algorithm we 
have adopted requires that the input data is positive, we 
have inverted the absorption-line profiles and set any nega- 
tive points to zero. For further details of the Roche tomog- 
raphy process see the review bv lDhillon fc WatsonI lEoOli) . 

Several Doppler imaging studies have used 

a tw o-temperature or filling-factor model (e.g. 

ICoUie r Cameron fc Unruhl |T994) during the image re- 
construction process. These studies work on the assumption 
that there will only be two-temperature components across 
the surface of the star, those of the immaculate photosphere 
and those of the cooler spots. Each image pixel is then ef- 
fectively made up of a combination of photosphere and spot 
intensity. The benefits of using a two-spot model are that 
it 'threshol ds' the maps, preven ting the growth of bright 
pixels (e.g. |H^ze^&^og^_^9^^jand preserves the total 
spot areas JCoUier Cameronl Il992h making quantitative 
analysis of the maps more easy. 

With Roche tomography we may well expect the sec- 
ondary star to display a much more complex topogra- 
phy than isolated stars. In particular, the surface tempera- 
ture across a CV secondary will vary greatly not just due 
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Figure 2. The LSD profiles and maximum entropy fits (see Section |3 for further details) for AE Aqr for observations starting on the 
night of 2001 August 9. The velocity of the secondary star has been removed using the parameters found in Section 1^1 and each profile 
is shifted vertically for clarity. The orbital phase for each exposure is indicated at the top right of each profile. 



to starspots but also as a result of irradiation, which is 
prominent in these systems even in single-line studies (e.g. 
i Davev fc Smitblll99a iDavev fc Smitbl 119961: IWatson et all 
l2003fl . For this reason, we cannot restrict our Roche tomo- 
grams to a two-temperature reconstruction. Thus Roche to- 
mograms present images of the distribution of line flux on 
the secondary star - effectively showing the contrast between 
the various features, blurred by the action of the regularisa- 
tion constraint. 



5 EPHEMERIS 

A new ephemeris for AE Aqr was determined by cross- 
correlation with a template star of spectral type K4V. We 
only considered the spectral region containing the strong Mg 
II triplet absorption for simplicity. The AE Aqr and K4V 
template spectra were first normalised by dividing by a con- 
stant, and the continuum was then subtracted off using a 
third order polynomial fit. This procedure ensures that the 



line strength is preserved along the spectral region of inter- 
est. 

The template spectrum was then artificially broadened 
to account for the rotational velocity (wsini) of the sec- 
ondary star. The projected rotational broadening of the 
secondary star was measured using an optimal-subtraction 
technique in which the template spectrum was broadened, 
multiplied by a constant and then subtracted from an 
orbitally-corrected AE Aqr spectrum. The process was re- 
peated in an iterative process, adjusting the broadening at 
each step, until the scatter in the residual spectrum was 
minimised. 

Through the above process a cross-correlation function 
(CCF) was calculated for each AE Aqr spectrum over the 
wavelength region covering the strong Mg ll triplet absorp- 
tion. The radial velocity curve in Fig. |S|was then obtained 
by fitting a sinusoid through the CCF peaks to obtain a new 
zero-point for the ephemeris of 

To = HJD 2452131.31345 ± 0.00007 (1) 
with the orbital period fixed at P = 0.41165553 d (from 
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Figure 3. The same as Fig. I^Jbut for the second night of WHT observations. The seeing deteriorated to ~2 arcsec towards the end of 
the night, hence the poorer signal-to-noise in the latter profiles. 



ICasares et al.l Il996l) . All subsequent analysis has been 
phased with respect to this new zero-point. This method of 
deriving the ephemeris is relatively insensitive to the use of 
an ill-matching template, or incorrect amount of broadening. 
Indeed, it is more likely to be dominated by the effects of ir- 
radiation, which is well known to introduce systematic errors 
in radial velocity studies if not properly acc ounted for (e.g. 
iDavev fc Srnitb[ll992l: IPavev fc Smithlll99fih . Although the 
steeper slope of the radial velocity curve around phase 0.5 
suggests the presence of irradiation which is confirmed in our 
Roche tomograms (albeit at a low level - see section [TJ, the 
ephemeris derived here provides a substantially improved 
image quality (both in terms of final reduced a^nd arte- 
facts in the reconstructed map) when compared against that 
achieved when using the zero -point of the ephemeris pub- 
lished bv lCasares et all l)l99(f) . 

We should note that the work in this section was carried 
out for the sole purpose of determining a revised ephemeris 
in order to improve the quality of the Roche tomograms in 
section Q From this perspective, no detailed attempt was 
made to determine the best-fitting spectral type and accu- 
rately measure the rotational broadening and other param- 



eters from this analysis - a more traditional study along 
these lines will be presented in a future paper. For com- 
pleteness, however, we obtained a secondary star radial ve- 
locity amplitude of Kr = 168.4 ± 0.2 km s~^ and a systemic 
velocity of 7 = -59.0 ± 1.0 km s"'^ (assuming the K4V tem- 
plate st ar HD24916A has a systemic velocity of -|-6.2±1.0 
km s~^. iKine et al.l 1200.^ 1. As discussed later in sectional 
these measurements will be biased due to, for example, sur- 
face inhomogeneities across the secondary star and the in- 
correct treatment of the absorption-line profile shape (which 
assumes that the star is spherical). Thus, the parameters de- 
rived from the radial velocity curve have not been used in 
the subsequent analysis presented in this paper. 



6 SYSTEM PARAMETERS 

As discussed in IWatson fc DhillonI lj200lfl . adopting incor- 
rect system parameters such as systemic velocity, component 
masses and orbital inclination when carrying out a Roche to- 
mography reconstruction results in spurious artefacts in the 
final map. These artefacts normally manifest themselves in 
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Figure 6. The measured secondary star radial-velocity curve 
witli tlie best-fitting sinusoidal fit for 2001 August 9 and 10. The 
data have been phase folded and then repeated. 

the form of bright and dark streaks and cause a decreased 
(more negative) entropy to resuh in the final map. fn other- 
words, adopting incorrect parameters causes more structure 
to be mapped onto the secondary star. By carrying out re- 
constructions for many pairs of component masses (iterating 
to the same on each occasion) we can construct an 'en- 
tropy landscape' (see Fig.|5J, where each grid element corre- 
sponds to the entropy value obtained in the reconstruction 
for a particular pair of component masses. One then adopts 
the set of parameters that gives the map containing the least 
structure (the map of maximum entropy). 

For the reconstructions carried out here we have as- 
sumed a root-square limb-darkening law of the form 

n^l)=I{l)[l-a{l-ll)-b{l-^)] (2) 

where /i = cos 7 (7 is the angle between the line of sight and 
the emergent flux), and 1(1) is the monochromatic specific 
intensity at the center of the stellar disc. The limb-darkening 
coefficients a = 0.647 and 6 = 0.213 were adopted, which 
is suitable for a star of T^f f — 4800K, log q = 4.5 at the 
central wavelength of our observations iClaretlll998f) . This 
form of the limb-darkening law has been chosen over the 
more usual linear dependency since it has been found that 
model atmospheres of cool stars show limb-d arkening laws 
which are largely non- linear fe.g. IClarelill998h . 

To determine the binary parameters of AE Aqr, we have 
constructed a sequence of entropy landscapes like that de- 
scribed above for a range of orbital inclinations i and sys- 
temic velocities 7. Then, for each i and 7, we picked the 
maximum entropy value in the corresponding entropy land- 
scape. From Fig. [7| this method yields a systemic velocity 
of 7 = -63 ± 1 km s~^, where we have adopted the 1 km 
s~^ error bar from the resolution of our grid search as a 
guide rather than a rigorous error e stimate. This i s con - 
sistent with the systemic velocity of ICasares et alJ (Il996l) 
who found 7 — —60.9 ± 2.4 km s~^, and is in excellent 
agreement with the va l ue of -63 ± 3 km s~^ found by 
I Welsh. Horne fc Gomerl 1^9^. The value obtained using 
Roche tomography is slightly more discrepant with the value 
of 7 = —59 ± 1 km s~^ obtained in section |^ from the ra- 
dial velocity curve. We note that the value of the systemic 
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Figure 7. Points show the maximum entropy value obtained in 
each entropy landscape as a function of systemic velocity, as- 
suming an orbital inclination of 66°. (The optimal value for the 
systemic velocity is, however, found to be independent of the as- 
sumed inclination). The solid curve shows a fourth order poly- 
nomial fit through these points with the intention to show more 
clearly where the peak in entropy occurs with respect to the as- 
sumed systemic velocity in the reconstructions. 
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Figure 8. Points show the maximum entropy value obtained in 
each entropy landscape for different inclinations, assuming a sys- 
temic velocity of 7 = -63 km s~ ^ . The solid curve shows a third or- 
der polynomial fit through these points. The entropy value peaks 
at an orbital inclination of i = 66° . 

velocity obtained from the entropy landscape technique is 
constant o ver a wide range of a ssumed orbital inclinations, 
as found in lWatson et al.l ll2003f) . 

Fig. |H| shows the maximum-entropy value as a func- 
tion of inclination obtained in entropy landscapes assuming 
the systemic velocity of 7 = —63 km s~^ derived above. 
This indicates a distinct increase in the amount of struc- 
ture present in the maps either side of the best fitting in- 
clination at i = 66°, which we have adopted as the optimal 
value for the orbital inclination of AE Aqr. It is comforting 
to see that this lies below the i < 70° limit inferred from 
the lack of eclipses : (.'lianaii. Middlcditch & Nelson 197Q), 
though it is still relatively h igh compared to previous es- 
timates. ICasares et alJ l)l996h determined i = 58° ± 6° by 
modelling the phase dependent variation of the rotational 
broadening in their intermediate resolution spectra, while 
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Figure 9. The entropy landscape for AE Aqr, assuming an or- 
bital inclination of i = 66° and a systemic velocity of 7 = —63 
km s~^. Both nights data were combined when carrying out the 
reconstructions. Dark regions indicate masses for which no ac- 
ceptable solution could be found. The cross marks the point of 
maximum entropy, corresponding to component masses of Mi = 
0.74 M0 and M2 = 0.50 Mq. 

I Welsh et al.l il995h constrained i to 55° ± 7° in their ab- 
sor ption line study of AE Aqr. 

[Sh^ib^ |l99^ noted, however, that the m ethods used 
bv lWelsh et alJ ^1995^ and lCasares et al] lll996ll to measure 
the orbital inclination of AE Aqr are prone to systematic er- 
rors and result in a shift to lower values for the inclination. 
This systematic error results from the incorrect treatment 
of the absorption-hne profile shap e which is often approxi- 
mated by using the Gray profile llGravlll992l) for spherical 
stars. Since the determination of the radial velocity ampli- 
tude and rotational broadening of the secondary star is also 
sensitive to the correct treatment of the line-profile - use 
of the Gray profile or failure to account for surface inhomo- 
geneities on the secondary star can lead to systematic errors 
in the derived mass ratio as well. Roche tomography, how- 
ever, predicts the correct broadening function for a Roche- 
lobe filling star and, in addition, takes into full account the 
line flux distribution across the star which can lead to sys- 
tematic errors in these more conventional studies. Further 
support for the higher inclination found in this work comes 
from a r ecent preliminary analysi s of high spectral- resolution 
data bv lEchevarn'a et al.l l|2005|). These authors have also 
found evidence for a high (~70°) inclination from studying 
the phase-dependent variation of the rotational broadening. 

The entropy landscape for AE Aqr assuming an orbital 
inclination of i = 66° and a systemic velocity of 7 = -63 
km s~^ is shown in Fig. 1^1 From this we derive a secondary 
star mass of M2 = 0.50 Mq and a primary mass of Mi — 
0.74 Mq. These are in good agreement with the most recent 
previous estimates in the literature, in cluding Mi = 0.79 ± 
0.16 M0 and M2 = 0.50 ± 0.10 Mq dCasares et alJll99et) . 
as weU as Mi = 0.89 ± 0.23 Mq and M2 = 0.57 ± 0.15 Mq 
JWelsh et al.iri995,l . 

In all cases we have carried out reconstructions to ~ 
0.8, though we have explored the effects of carrying out the 
reconstructions to different values. For reconstructions to 
lower values the entropy landscapes break up into noise. 



whereas when iterating to higher x^ values the entropy land- 
scapes favour higher inclinations (which are unphysical given 
the lack of eclipses) - eventually losing all sensitivity to the 
assumed inclination. In all cases, however, the derived com- 
ponent masses for any given orbital inclination all agree to 
better than O.O2M0 for all reasonable selections of target x^- 
We are confident that we are neither over-fitting nor under- 
fitting the data. The maps rapidly break into noise for x^ 
values less than 0.8, leading to noisy entropy landscapes. 
Maps fit to greater x^ values become 'washed out' with fea- 
tures becoming less distinct as more pixels are assigned the 
default value in the absence of proper data constraints. 

It should also be noted that we have not assigned er- 
rors to our system parameter estimates. Technically, this 
could be done using a Monte-Carlo style te chnique and 
applying the bootstrap resampling method jEfronl Il979l : 
^fron & Tibshirani 1993) to generate synthetic datasets 
drawn from the same parent population as the observed 
dataset (see also IWatson fc DhillonI [2OO j) . Unfortunately, 
this would mean repeating the analysis carried out in this 
paper for several hundred bootstrapped datasets which 
would require several months of computational time. It is, 
therefore, unfeasible for us to assign strict errors bars to our 
derived binary parameters. Given the quality of the data 
(high time and spectral-resolution) and the full treatment 
of the data by Roche tomography (which takes into account 
the roche-lobe shape and any surface inhomogeneities) , we 
believe that the system parameters derived in this paper for 
AE Aqr are the most accurate to date. 



7 SURFACE MAPS 

Using the binary parameters derived in Section |S| and com- 
bining both nights datasets, we have constructed a Roche 
tomogram of the secondary in AE Aqr, and this is displayed 
in Fig. 1101 The corresponding fit to the datasets obtained 
on both nights are displayed in Fig. 1111 Dark spots are ev- 
ident in the Roche tomogram. The most noticeable feature 
is the large starspot feature which lies near (but not centred 
on) the polar regions at a latitude of ~65° and is skewed 
towards the trailing hemisphere ((/)=0.25) of the star. Such 
high-latitude and polar spots are common in Doppler images 
of other rapidly-rotating stars and are discussed in more de- 
tail in Section |S] 

The other prominent dark region in the tomogram is 
best seen at 0=0.5. This lies at the Li point on the sec- 
ondary star, and could possibly be due to irradiation from 
the accretion regions which would also appear dark in the 
tomograms since the weak metal lines become ionised. Al- 
though irradiation of the donor star in AE Aqr has not been 
seen previously, it's location at the Li point is consistent 
with that seen in other CVs (e.g. Watson et al. 2003). Fur- 
thermore, in high spectral-resolution observations taken at 
the Magellan 6.5-m telescope during July 2004 (to be pre- 
sented in a later paper), the preliminary LSD profiles around 
0=0.5 are almost identical to those presented here. This is 
indicative of a fixed feature around the Li point which is 
more likely due to the presence of irradiation rather than 
the existence of a long-lived (several year) spot feature. On 
the other-hand, we currently do not have a knowledge of the 
spot life-times for this class of star and the modelling of the 
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mass transfer history of AM Her bv lHessman et all (|200(J) 
suggests that the spottedness of the Li point on the donor 
stars in CVs may be untypical compared to isolated stars. 
Despite these caveats, however, irradiation is still our pre- 
ferred interpretation of the feature located at the Li point. 

In addition to these prominent features, a number of 
other smaller spots can be seen at lower latitudes. Un- 
fortunately many of these lower latitude spots are heav- 
ily smeared across the equator since the latitudinal reso- 
lution of the imaging technique decreases towards the equa- 
tor (this is because the path of a starspot trail in the time 
series shows the least variation as a function of latitude 
near the equatorial regions). This problem is further com- 
pounded by the lower signal-to-noise of our LSD profiles 
relative to many other Doppler imaging studies of isolated 
stars (a result of the faintness of AE Aqr), and the tar- 
get's relatively high orbital inclination, which results in a 
degener acy between which hem isphere a feature is located 
on fe.g. lWatson fc Dhillonll200ir) . Therefore, it is difhcuh to 
tell whether the lower latitude spots are located within a 
particular latitude band or not. 

To check whether or not the features seen in the Roche 
tomogram are real, or spurious artefacts due to noise, we 
have carried out two further Roche tomography construc- 
tions using only the odd-numbered spectra for one and the 
even-numbered spectra for the other (see Fig. I13II . Despite 
the reduction in phase-sampling, both of these Roche to- 
mograms show the same features as seen in Fig. EOl and 
therefore we are confident that these features are real. Fur- 
thermore, when we construct separate maps for each nights 
data (not shown here for reasons outlined below), we again 
find that that features are reconstructed at similar locations 
in the maps where there is sufficient phase overlap in the 
observations. This gives us further confidence in the reality 
of the reconstructed features. 

Although the maps reconstructed on individual nights 
show similar features, these features do differ in detail. Un- 
fortunately, this is most likely a result of the relatively poor 
signal-to-noise of our observations combined with the in- 
troduction of phase gaps when reconstructing each dataset 
separately - further compounded by the deteriorating seeing 
conditions on night 2. We feel, therefore, that it would be 
unwise to speculate as to whether any of the detailed dif- 
ferences between the individual nights reconstructions are 
due to the evolution of spots other than to mention that 
there does not appear to be any significant movement of 
spot features. Future higher signal-to-noise observations will 
be better placed to resolve spot evolution. 

In order to make a more quantitative estimate of the 
spot parameters on AE Aqr, we have examined the inten- 
sity distribution of the pixels in the Roche tomogram. By 
looking for a bimodal distribution in pixel intensities, an es- 
timate of the spot properties can be obtained by labelling 
low-intensity pixels as spotted regions, and high-intensity 
pixels as immaculate photosphere. First, we discarded all 
pixels on the southern hemisphere as this hemisphere is least 
visible and as such will mainly contain pixels assigned with 
the default intensity which, in this case, will be the average 
intensity of the map. Inclusion of these pixels would 'wash- 
out' any underlying bimodal distribution in pixel intensity. 
The corresponding histogram of pixel values is shown in 
Fig. 1121 and, as expected, this shows a group of pixels of 



lower intensity which have been labelled as starspots (see 
caption of Fig. ll2l for details). 

Although far from an ideal means of classifying what 
regions are spotted and what regions contain immaculate 
photosphere, this is the best classification scheme given our 
inability to apply a two-temperature model as discussed in 
Section 0] Based on this scheme. Fig. 1141 shows the distri- 
bution of starspots as a function of latitude scaled by either 
the total surface area of the star, or by the surface area of 
the latitude strip over which the spot coverage was calcu- 
lated. This confirms our visual inspection of the tomogram 
and that the density of spots increases towards the pole, 
with a grouping of spots centred around a latitude of -1-65° 
corresponding to the large spot seen near the north pole. 

Inspection of Fig. 1141 also shows an apparent depletion 
of starspots near -1-40° , with an increase in starspot coverage 
again towards lower latitudes, with a further increase in spot 
coverage below -1-10°. This latter increase in spots below a 
latitude of ~ 10° is probably due in part to the presence 
of irradiation near the Li point. Nonetheless, the irradiated 
region is not large enough to fully explain the increase in 
low-latitude pixels being flagged as spots. Therefore, this 
indicates the presence of a lower latitude site of magnetic 
activity on AE Aqr which future higher signal-to-noise ob- 
servations will help resolve. 

In total, we estimate that 18 per cent of the north- 
ern hemisphere on AE Aqr is spotted. We stress that this 
number should be regarded with caution given our scheme 
for classifying pixels containing spots, but that such a fig- 
ure for the spot filling-factor is not at all unreasonable 
and, if anything, is likely to be an under-estimate of the 
tru e spot filling f a ctor. TiO studies of rapidly rotating stars 
by lO'Neal et al.l il998l) have revealed spot filling factors 
as large as 50 per cent in the case of II Peg. Generally, 
such studies result in much larger filling factors than are 
deduced from Doppler images which typically return spot- 
fiUin g factors from ~4 p er cent (e.g. Speedy Mic, 4.5 per 
cent. lBarnes et al. 2001811 to ~ 12 per cent (e.g. He 699, 11.9 
per cent. lBarn£^^^2001hll . 

A prime example of the underestimation of the global 
sp ot coverage by Dopp ler imaging methods was highlighted 
bv lBarnes et alJ l)2004h . whose Doppler images of the pri- 
mary star in the contact bi nary AE Phe return ed a spot 
filling factor of ~5 per cent. lBarnes et al] l|200^, however, 
estimated that the actual spot coverage could be 27 per cent 
in order to lower the apparent temperature of the primary 
component to the level seen in their observations. This dis- 
crepancy between Doppler imaging and, for example, TiO 
studies is probably caused by the failure of Doppler imag- 
ing to resolve small starspots, which will certainly be the 
case for the lower signal-to-noise dataset presented in this 
paper. The poorer latitudinal resolution near the equator 
will also cause low-latitude spots to be heavily smeared over 
a larger surface area, reducing their contrast with the im- 
maculate photosphere and increasing the probability of the 
non-detection of s tarspots at these lower latitudes. Finally, 
IWebb et alJ ll2002ll found a spot filling factor of 22 per cent 
from a TiO study of the CV SS Cyg, similar to our estimate 
for AE Aqr. 
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Figure 11. Roche tomography fits to the LSD trailed spectra showing the observed trails, computed trails from the Roche tomography 
reconstruction and the O-C residuals. The greyscale wedge (right) indicates the magnitude of the residuals as a fraction of the average line 
profile depth. The largest residuals result from spikes in the deconvolved continuum (which do not adversely affect the Roche tomography 
fits to the profiles themselves - see Figs|5]&|2J. The residuals are greatest for the second nights observations due to the deterioration in 
the seeing at the end of the night. Top: observations from 2001 August 9. Bottom: observations from 2001 August 10. 
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Figure 13. As for Fig. 1101 but reconstructions for odd-numbered spectra (left) and even-numbered spectra (right). 
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Figure 14. Plots showing the spot coverage as a function of latitude. Left: spot coverage as a function of latitude expressed in terms 
of the surface area at that latitude. Right: spot coverage as a function of latitude normalised by the total surface area of the northern 
hemisphere. The curves appear noisy in comparison to similar plots elsewhere in the literature as the elements that make up our stellar 
grid are not aligned along constant latitude strips in order to ensure that the inner-Lagrangian point is tiled adequately. Therefore a 
certain degree of interpolation in latitude between adjacent tiles is required to calculate the contribution of each tile to that latitude. 



8 DISCUSSION 

The secondary in the cataclysmic variable binary AE Aqr 
is a highly active star as evinced by the heavily spotted 
Roche tomogram presented in this work. This is unsurpris- 
ing, given the rapid (0.41-d) rotation period of AE Aqr. In 
keeping with the abundant observations o f high-latitude or 
polar spots on rapidly rotating stars (see ISchriiver fc Titlel 
I2OOII and references therein), AE Aqr also exhibits a large 
spot at high-latitude. This spot is, however, not truly polar 
as seen in many Doppler images as it is centred at a lat- 
itude of around +65°. This is still in striking contrast to 
our Sun, where s tarspots seldom appea r beyond ±30° from 
the equator (e.g. iPulkkinen et alJll999l) . After much debate 
about the re ality of the polar spots revealed by Doppler 
imagmg (see 'Hatzes et aTTll99el. lUnruh fc Collier CameronI 
[1997 , Bruls et al. 1998), the reality of t hese features now ap- 
pear t o be established within the field. ISchiissler fc Solankl 
l|l992() have suggested that such high latitude spots are 
the cause of the strong Coriolis forces in these rapidly ro- 
tating stars which drive the magnetic flux tubes towards 



the polar regions. Alternatively, such spots can form at 
lower latitudes and then migrate polewards - such a pole- 
ward migration of mid-latitude spots has been reported 
for the RS CVn binary HR1099 (V711 Tauri - Vogt et al] 
ll999llStrassmeier fc Bartujl200d) . Future observations will 
be needed in order to establish whether the starspots on AE 
Aqr also undergo a poleward motion. 

In addition, AE Aqr also appear s to show starspots a t 
latitudes below ~ 35°. As described in lCranzer et al.li200Cl) . 
such bimodal spot distributions (the simultaneous appear- 
ance of spots at high and intermediate/low latitudes) have 
also been reported for the zero-age main-seq uence (ZAMS) 
rapid rotator AB Dor (e.g. lUnruh e t aL* 1995), the Ploiades- 
age ma i n-seq uence star LQ Hva fclice fc Strassmeier 1993 . 
iDonati Il999l) , the a Per G-dwarves He669 and He520 
l^anMset^J; 'l99 ^, and the young G-dwarf EK Dra 
(Strassmeier fc Ricel ll99'3) amongst others. 

In the ir paper studying the distribution of starspots on 
cool stars, ICranzer et all ll2000t) found that low -mass stars 
(~ 0.4Mq) with rotation rates more than 4 times solar 
showed spot emergence at latitudes between ~ 60 — 70°. 
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Figure 10. Roche tomogram of AE Aqr. Dark greyscales depict 
regions of suppressed absorption and are indicative of starspots 
or irradiated zones. The system is plotted as the observer would 
see it at an orbital inclination of 66° , except for the central panel 
which shows the system as viewed from above the North pole. 
The orbital phase (with respect to the ephemeris of Section |SJ is 
indicated above each panel. For clarity, the Roche tomograms are 
shown excluding limb-darkening. 
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Figure 12. Histogram of the pixel intensities from the Roche 
tomogram of AE Aqr after pixels on the southern hemisphere 
(latitude < 0°) were discarded. The brightest pixel in the map is 
assigned an intensity of 100 and other pixel intensities are scaled 
linearly against this. For the purposes of estimating the global 
spot properties of AE Aqr, the immaculate photosphere has been 
defined as pixels with intensities of 74 or greater, and the spotted 
photosphere those regions where pixel intensities are less than 74. 



It is perhaps worth noting that this is comparable with AE 
Aqr which, at O.5M0 , shows spots emerging at ~ 65° . Unfor- 
tunately, since AE Aqr has already evolved some way off the 
ZA MS, compariso n to th e results from modelling young stars 
by iGranzer et alj ||200(J) may make any conclusion drawn 
from such a comparison somewhat uncertain. 

It is also interesting to compare the images of AE Aqr to 
those of the late pre-main sequence o r early ZAMS isolated 
star Speedy Mic jBarnes et al. 2001all . Excluding their ages 
and the effect of tidal forces, Speedy Mic (of spectral-type 
K3V and rotation period 0.38-d) otherwise has similar stel- 
lar parameters to AE Aqr {Prot = 0.41-d, K3-5V/IV). Like 



AE Aqr, Speedy Mic does not show any truly polar spot 
but does show a concentration of starspots around 60-70° 
at a similar latitude to the large spot group observed on AE 
Aqr at 65°. Speedy Mic also exhibits a bimodal distribu- 
tion of starspots, with a similar concentration of spots at 
lower latitudes to that observed for AE Aqr. Furthermore, 
Barnes et al. (2001a) also describe a relative paucity of spots 
at a latitude of 30-40°, the upper limit of which coincides 
with the depletion of spots around 40° seen on AE Aqr. 
Such depletion of spots at interm ediate latitudes is also ob- 
served in the K2 dwarf LQ Hya iRice fc StrassmeierlllQQ^ . 
Donatill999t) and the K5 dwarf BP 4- 55° 4409 (LO Peg - 
Lister. Collier Cameron fc Bartujll999^ . More observations 
of late-type stars are required in order to determine whether 
this feature is common to all K-type and later stars. 

In conclusion, we have unambiguously imaged, for the 
first time, large starspots on the secondary star in AE Aqr, 
the distribution of which is similar to other rapidly rotat- 
ing low-mass stars. We estimate a spot coverage of approx- 
imately 18 per cent, indicating a high-degree of magnetic 
activity which is no doubt due to an efficient dynamo driven 
by rapid rotation. Such a high-level of activity on AE Aqr 
provides evidence that magnetic-braking in long-period CVs 
could play a significant role in their evolution towards the 
period gap. 



9 FUTURE OPPORTUNITIES 

Given that this is the first time that starspots have been 
imaged on the secondary stars in CVs, we thought it prudent 
to discuss the future avenues that such studies could explore. 



9.1 Spot properties on CV secondaries 

As briefly mentioned in Section a knowledge of the dis- 
tribution of starspots on Roche-lobe filling stars can pro- 
vide tests of our under standing of stellar dynamos. For ex- 
ample, the models of Holzwarth fc Sch iissler (2003) show 
that the tidal effects in a binary due to the presence of a 
close companion can lead to the formation of clusters of flux 
tubes (and therefore starspots) at preferred longitudes. Fur- 
thermore, these effects are predicted to become more pro- 
nounced for shorter period systems where the tidal forces 
are strongest. Since CVs contain some of the most rapidly 
rotating stars known, they provide excellent candidates in 
which to study the effects of tidal forces on magnetic flux 
tube emergence. In addition, since the magnitude of the tidal 
effects will change with changing binary parameters, so will 
the predicted starspot distributions. Therefore it is highly 
desirable to obtain a set of high-resolution maps of a num- 
ber of different CVs as each individual binary system stud- 
ied should, in principle, provide an independent test of such 
models. 

Studies should also be carried out to investigate how the 
starspot coverage varies with time in an attempt to deduce 
whether these stars display magnetic activity cycles similar 
to the 11-year solar cycle. An understanding of these cy- 
cles is crucial to the success of any stellar dynamo theory. 
As outlined in Section magnetic activity cycles have also 
been invoked to explain quasi-cyclical variations in orbital 
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periods, quiescent magnitudes and outburst intervals, dura- 
tions an d shapes, with cycles on t imescales as short as ^ 
6 years l|Ak. Ozkan fc Matteilliool . No direct evidence for 
magnetic activity cycles has so far been found for these sys- 
tems and therefore long-term monitoring of a selection of 
CVs is actively encouraged. To this end we have initiated a 
long-term monitoring campaign of AE Aqr, the first results 
of which are to be published in a future paper. 

9.2 Differential rotation of tidally distorted stars 

Differential rotation is thought to play a crucial role in am- 
plifying and transforming the initial poloidal magnetic field 
into toroidal field through dynamo processes (the so-called SI 
effect). Thus measurements of differential rotation rates aim 
to evaluate the influence of stellar parameters on the surface 
shear rates and hence test stellar dynamo theories. Theoret- 
ical wo rk regarding d i fferen tial rotation in tidally distorted 
stars bv lScharlemannI lll982ll suggests that it should be weak- 
ened, but not entirely suppressed, in these objects. Some 
two decades later, however, observational evidence for the 
strength of differential rotation in tidally distorted stars is 
still confused. Petit ct al. (2004) find evidence for weak dif- 
ferential rotation in the RS C Vn system HR 1099 . Images 
of the contact binary AE Phe iBarnes et aljl20o3) . on the 
other-hand, shows rapid spot evolution and motions. Evi- 
dently more studies of tidally distorted systems are required. 

By acquiring two images of a CV separated by several 
days, one can cross-correlate strips of constant latitude, al- 
lowing the sense and magnitude of the differential rotation 
to be measured. Furthermo re, in the event that surface shear 
exists. IScharlemannI (Il982l) note that, on average, the stellar 
envelope should still be forced to co-rotate. Roche tomog- 
raphy studies could also be used to deduce the co-rotation 
latitude at which the envelope is tidally-locked to the orbital 
period. 

Although fainter than RS CVn systems, CVs still of- 
fer an excellent laboratory to study differential rotation in 
tidally distorted stars. For instance, RS CVn's have orbital 
periods that typically span several days, meaning that long 
observations covering several nights are required to estab- 
lish a complete picture of the system. During this time in- 
dividual starspots may have appeared, evolved into differ- 
ent shapes or disappeared entirely. Therefore a single im- 
age of an RS CVn system can often not be regarded as a 
single 'snapshot' of the spot morphology - making the de- 
tection and accurate determination of differential rotation 
rates more complicated. CVs, however, have shorter orbital 
periods (typically < 10 hours) and therefore a single Roche 
tomogram of a CV should be obtainable in one night for the 
majority of cases - making differential rotation studies more 
accessible. 



9.3 Impact of activity on accretion dynamics 

Magnetic activity on the secondary stars in CVs is also 
thought to play a role in the accretion dynamics of in- 
teract ing binaries over short timescales. iLivio fc Pringld 
il994) have suggested that the low states observed in 
some interacting binaries can be explained by the passage 
of starspots across the inner-Lagrangian (Li) point (see 



also lKing fc Cannizzdll993) . This temporarily quenches the 
flow of material through the Li point, thereb y reducing 
the accretion rate. Indeed. iHess man et alJ ll2000ll were able 
to model the long-term brightness variations of the cata- 
clysmic variable AM Her as modulations caused by starspots 
and conclude that either these systems have a large spot 
covering fraction (~ 0.5) or that the Li point i s unusu- 
ally spotty. More recently, Ide Martino et aO ll2002l) also at- 
tributed the rapid variation of accretion in Am Her to 
starspots. Such mass transfer variations can also have a 
marke d effect on the respo nse of the accretion disc in dwarf 
novae dSchreiber et al.llioO O'') . 

There are a number of ways that these ideas can be 
tested with Roche tomography. First, the observation of a 
starspot crossing the Li point in combination with a simul- 
taneous drop in system brightness (high-resolution Roche 
tomography datasets should normally be accompanied by 
simultaneous photometry) would provide evidence support- 
ing these models. Although such an observation may seem 
fortuitous, given the large spot covera ges seen on rapidly 
rotating stars (as much as 50 per cent, lO'Neal et alJll99dl 
and, in particular, given the high spot coverages seen in 
the Roche tomograms of AE Aqr, there is a high proba- 
bility of witnessing such an event. The main problem with 
this approach would be assessing whether a feature at the 
Li point was truly a spot or caused by irradiation. One 
suitable target is t he bright Z Cam type dwarf-nova V426 
Oph (iHellier et ahi riggO) which exhibits drops from stand- 
still which may be due to a decrease in mass transfer. (In 
recent high-resolution Magellan observations starspot fea- 
tures were found in the deconvolved line profiles of V426 
Oph and a Roche tomogram of V426 Oph will be published 
in a later paper). 

Second, even if an event is not seen directly, given a de- 
tailed knowledge of the starspot distributions and the sense 
and magnitude of their motions of the starspots due to dif- 
ferential rotation (Section 19.211 . it should be possible to pre- 
dict at what times a starspot will migrate across the Li 
point. Again, an observed dimming at the predicted times 
would pro vide strong eviden c e in fa vour of such models. Fur- 
thermore, jKin^^^Camizz^ il998l) note that the most rapid 
transitions to low states occur over the course of a day. The 
leng th scale of the mass overflow region near the Li point is 
fe.g. lPringljll985D 

w~^~3x lO^Phr cm, (3) 
27r 

where Cs is the sound speed of the gas and Phr is the orbital 
period P expressed in hours. In order to cause a transition 
to a low-state over one day, a spot boundary would have 
to move across the Li point at a velocity of ~ 3 x 10'^ Phr 
cm s~^. By using starspots as surface flow tracers it will be 
possible to determine whether the differential rotation rates 
in CVs are commensurate with these velocities. 

Finally, in the model for AM Her proposed by 
iHessman et al.l (l2000l) . they infer a power- law for the dis- 
tribution of spot sizes. Using Roche tomography, one can 
determine both the numbers and sizes of the starspots on 
CV secondary stars, allowing the measurement of the power- 
law index for the uppe r end of the spot size distribution (e.g. 
ICollier Cameronll200lD . This would provide a powerful con- 
sistency test for such models. 
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In conclusion, not only are Roche tomography studies 
important for our further understanding of the diverse be- 
haviour of CVs but they also promise to provide further 
insights into stellar dynamos. In order to provide the best 
constraints on spot sizes and distributions, and to extend 
surveys to more (and almost certainly fainter) CVs, observa- 
tions using echelle spectrographs on 6 to 8-m class telescopes 
will be required. 
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